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Abstract 

a! 

Ph! We use thermodynamic Bethe ansatz to study nonrelativistic scattering theory of 



low energy excitations of ID Hubbard model, using the S'-matrices proposed by Efiler 
and Korepin. This model can be described by two types of excitation states, holons 
and spinons, as asymptotic states. In the attractive regime, the spinon is massive while 
the holon is massless. The situation is reversed with a repulsive coupling. We derive 
that the central charge of the Hubbard model in the IR limit is c = 1 while it vanishes 
in the UV limit. The contribution is due to the massless degree of freedom, i.e. the 
holon for the attractive regime, and the massive mode decouples completely. This result 
is consistent with various known results based on lattice Bethe ansatz computations. 
Our results make it possible to use the S'-matrices of the excitations to compute more 
interesting quantities like correlation functions. 
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There has been considerable interest in the interplay of integrable quantum field 
theory and statistical mechanics 0. In particular, a lot of progress in this relationship has 
been made in two-dimensional models. One of the most useful methods in these models 
is the factorizable S-matrix theory [0]. In 1 + 1-dimensional integrable field theories 
where infinite number of conservation laws exist, the S- matrices are facorizable into two- 
body elastic S-matrices, which satisfy the Yang-Baxter equation. With known particle 
spectrum and additional symmetries, one can determine the S'-matrices exactly^] by 
solving the Yang-Baxter equation. Another method is to diagonalize the Hamiltonian 
using Bethe ansatz and find physical particle states and their S-Hiatrices||. In addition 
to their importance as physical amplitudes of scattering between asymptotic particle 
states, exact S'-matrices can give other interesting quantities like the central charges 
of underlying conformal field theories (CFT's) from the finite size effects, conformal 
dimensions of the operators, and even the correlation functions. 

Our motivation is to establish the S-matrix approach to study the non-relativistic 
lattice models like the ID Hubbard modeljS|. Although the Bethe ansatz method is 
quite useful in finding eigenvalue spectra, excitation states, and their thermodynamic 
properties, it is not so useful in finding other important quantities, in particular, cor- 
relation functions. Recently, Efiler and Korepin derived the S-matrices for low lying 
excitations of the ID Hubbard models . What concerns us here, as a first step towards 
the complete S-matrix bootstrap of the low lying excitations of the Hubbard model, is 
how to confirm the validity of these matrices. For this purpose, we employ thermody- 
namic Bethe ansatz (TBA) for 2D QFT which is now a standard tool@ to check the 
S-matrices. In the original formulation of factorizable scattering theory, one is inter- 
ested in the scattering of relativistic excitations. However, in interacting 1-D quantum 
systems, such as the Hubbard model, there are in general several low-energy excitations, 
with complicated dispersion relations and with different Fermi velocities. Also, the mo- 
tions of low energy excitations are far from relativistic. Thus, what we should do is to 
extend this TBA method to non-relativistic theories to compute the central charges in 
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the IR and UV limits and compare with results of finite size effects 0. 



Among the lattice models of the strongly correlated electron systems in low dimen- 
sions, it is believed that the two-dimensional Hubbard model provides some of the prop- 
erties of high-T c superconductivity 0. Furthermore, a strong quantum fluctuation in 
low dimensions suggests common features in 2D and ID Hubbard models. Fortunately 
ID Hubbard model can be exactly diagonalized via the Bethe-ansatz technique ||. Its 
thermodynamic properties, such as the susceptibility, the magnetization, and the low- 
temperature specific heat, for the repulsive (U > 0) and the attractive (U < 0) on-site 
interactions have been studied in the literature ||. It is noticeable that the excitation 
spectrum is described by the spin and charge excitations, i.e. spinons and holons, and 
the spin (charge) excitation possesses a gap in the attractive (repulsive) Hubbard model 
as long as the on-site interaction U exists. The low-energy charge (spin) excitations 
for the attractive (repulsive) case are proportional to those for the antiferromagnetic 
Heisenberg chain irrespective of the strength of the interaction and of the electron fill- 
ing. Woynarovich and Eckle|7| have analyzed the finite-size corrections in the Hubbard 
model for the repulsive and half-filled case, and their analysis yields the central charge 
of the Virasoro algebra c = 1 which is contributed by the spin excitations, while the 
contribution of the charge excitations is negligible only if the on-site repulsion U is not 
so weak. 

Let us now examine the dispersion relations of the excitations in the attractive case. 
Note that the results for the repulsive case can be obtained from the attractive case by 
interchanging the roles of holons and spinons. The holon energy in terms of rapidity A 
is given by 



In the above U is the coupling constant of the on-site interaction of the ID Hubbard 
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(1) 



whereas, momentum is given by 




model. Notice that the holon is massless and E c ~ v c P c , as A — > oo, where t> c , the Fermi 
velocity, is given by v c = E' c {oo) / P' c {oo) . 

Spinon has the following dispersion relations, for energy, 

00 duj J\{uj) cos(c<j sin k) 



r 

EJk) = 2\U\ -2 cos k + 2 / 

Jo 

and momentum, 
P s (k) = k 



uj cosh(ci;?7) 



exp(—\uU\), 



duo Jn(oj) sinfcj sin k) , . n n 
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We see that the spinon stays massive for finite U for all values of the rapidity. 
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Recently, using these dispersion relations, Efiler and Korepin derived the scattering 
matrices of Hubbard model from the Bethe ansatz solution, by generalizing a method 
of extracting S'-matrix from the asymptotics of the wave functions of the scattering 
state H]. These ^-matrices of the excitations on the lattice are well-defined as far as 
the wave packets of two excitations are well-separated. The complete scattering matrix 
is 16 x 16 dimensional and is in a block diagonal form consisting of 4 blocks. Each of 
the blocks describe scatterings of holon-holon, spinon-holon, holon-spinon, and spinon- 
spinon, respectively, as follows: 



S 



SrJu) 















S.Jw) 



S cs (w) 



S ss (v) 



(5) 



where the holon-holon scattering amplitude is 

r(±±^)r(i-f) 



s ~ u r (i^)r(i + f) 
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A - A' 
2\U\ 



(6) 
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with T being the gamma function. In the above I and P are the 4x4 identity and 
permutation matrices respectively. This S'-matrix has the same form as that of the spin 
| Heisenberg antiferromagnet and of the 577(2) i WZNW model|ll[]. The spinon-spinon 
scattering amplitude is 

c ( v ) - r (^) T ( 1 + ^) ( v i i F \ sinfc-sinfc' 

b " {v) - r (m)r(i-*) \v-i v-i J' 2\U\ ' {V 



2 J \ 2 

and can be obtained from S cc by setting u — > — 1>. The spinon-holon scattering amplitude 
is, 

.l+zexp(vrw) A -sin A; 

5 SC W = ; ? -I, W = . 8 

1 — iexp{nw) 2\U\ 
We have the same form for holon-spinon scattering amplitude S cs . Notice that S sc and 

S cs approach constant values as w — > oo. 

Let us now formulate the TBA for non-relativistic scattering. The TBA computes 
the Casimir energy of a theory on a circle of length R with S'-matrices and particle 
spectrum as input data0. With a temperature T — 1/R the configuration of minimizing 
free energy gives the ground state energy of the system, which is again related to the 
central charges of the underlying UV CFT by 

E smund (R)~-J2^ (9) 

i 1 

where V{ are Fermi velocities for excitations in the system and Cj are the corresponding 
effective central charges. 

Consider N particles in a box of length L with periodic boundary condition(PBC). 
Moving a k-th particle of type a with energy E a (8k) and momentum P a (0k) all the way 
by exchanging with other particles and coming back to the original configuration using 
PBC, we get 

JV 

e -iLP a {e k ) -q Saai (9 k -9i) = h (10) 

where the index Oj specifies species of the i-th particle. In general the product of S'- 
matrices is a large size matrix called transfer matrix and one should diagonalize this 



by some techniques. However, we consider a diagonal scattering theory first because 
the non-diagonal case can be understood as a slight modification. Taking logarithms on 
both sides of Eq.(lO) gives 

N , 

- LP a (6 k ) + £ - In S aai {9 k - 6i) = 2irn k , (11) 
i=\ 1 

with an arbitrary integer n k . In the thermodynamic limit, N, L oo, one can express 
Eq.(ll) in terms of density of states like 

2n Pa (6) = -lp'm + W de' P \{e')<t> ab {e - e'), (12) 

b J 

where p a {9) and p l a {9) are the densities of allowed and occupied states, respectively, and 
<f) a b is the logarithmic derivatives of S'-matrices S a b- In terms of 'pseudo-energies' e a 
defined by e~ ea = p\j (p a — pi), one can express the ground state energy by 

E gTomid (R) = ~ E /_~ ^RK(0)LeM, (is) 

where L e (9) = ln[l + e _e ^]. e a is determined by the minimizing condition of the free 
energy which is the following set of nonlinear equations: 

e a {9) = RE a (9) -Y^^ab* L £b (9), (14) 

b 

where * denotes rapidity convolution, f * g(9) = J^° 00 d9' f(9 — 9')g(9'). As we mentioned 
above, the sum in Eq.(13) for non-diagonal theories has to been taken with care. Di- 
agonalization of transfer matrix brings in additional 'massless' particles, which do not 
contribute to the central charge directly in the sum, due to masslessness, but which nev- 
ertheless affect the massive particle distributions. Additional care for non-relativistic 
scattering is that the participating particles can have different Fermi velocities, unlike 
in the relativistic case where all have the light speed as Fermi velocity. 

For the non-diagonal theories, the product of S'-matrices in Eq.(lO) is replaced by 
the eigenvalues of the transfer matrix. For ID Hubbard model the nondiagonal matrices 
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S cc and S ss are six vertex model type: 



S a 



a a 



b a c a 



c a b a 



V 



a r 



a = c, s. 



(15) 



The eigenvalues of the transfer matrices and the associated constraint equations can be 
derived by algebraic Bethe Ansatz method to be 

N M n / n\ N M 

a a {yr - t>) 
\b a (y r -6) ; 



K(0) = \{a a (6-6 l )Xl 



and 



i=l 



N 

n 



QalP ~ Vr) 
1 &a - yr) ' 



(16) 



6 a (2/fc - tt b a(yr - Vk) a a (y k - y r ) 



-1. 



. . ,i ■■ (17) 

^ a«(yfe - 0<) r= i a a (y r - y fc ) 6 Q (y fc - y r ) 

From the S'-matrices Eqs.(6) and (7), one can read off the corresponding matrix elements 

to evaluate the explicit eigenvalues. The holon and the spinon sectors are coupled by 

diagonal scattering matrix S sc . 



Hence we have two kinds of the periodic conditions for the Bethe wave functions of 
holons and spinons in terms of A c , A s , and S cs and two sets of the constraint equations. 
For simplicity, we concentrate on holon sector first and will generalize the argument to 
spinons. From Eqs.(6) and (17) the constraint equation for holon sector becomes 



7 Pk - ®i 



n 



i\U\ *^ p k -/3 r + 2i\U\ _ 1 



(18) 



t = \ Pk~0i + i\U\ r t\ Pk-Pr- 2i\U\ 

where we have introduced the shifted rapidities Pi — yi — i\U\ to have the unitary form. 
It is well known, from the analogy in the antiferromagnetic Heisenberg chain, that in the 
thermodynamic limit, iV — > oo, the general solutions of these equations are the strings 
consisting of n-pseudoparticles of roots = /3° + i\U\(n + 1 — 2j), where is real, 
j = 1, • • • , n, and n = 1, • • • , oo. Such a string is a bound state of n-pseudoparticles 



and can be interpreted as a fictitious massless particle of real rapidity /3°. Since the 
length of the strings could be infinitely long, there are infinite number of such massless 
particles. Similarly the constraint equations for spinons can be understood in the context 
of another kind of pseudoparticles designated by the rapidities a's which also form the 
string-solutions in the thermodynamic limit. Therefore what we have is a diagonal 
scattering theory of holons, spinons and infinite number of massless particles associated 
with them. The scattering amplitudes of n-th massless particle with holon is 

9-in\U\ 



Sn(0) n .,,, r 
d + in\U | 

and the scattering amplitudes between the massless particles are 
9 + i\n — m\\U\ 



(19) 



S n m{9) 



9 — i\n — m\ \U\ 
+ i(|n-m| + 2)|l7| 9 + i{n + m - 2)\U\ 



9 + i(n + m)\U\ 
9 -i(n + m)\U\ 



(20) 



9 - i(\n - m\ + 2)\U\ 9 - i(n + m - 2)\U\ 

For the spinons, the corresponding scattering amplitudes are obtained by the replace- 
ment 9 — > —9. 



The minimizing condition of the free energy using the above mentioned S'-matrices 
leads to an infinite set of non-linearly coupled equations. It is a standard procedure to 
use Fourier transformation on the TBA equations and to simplify the equations in terms 



of a unified kernel if = (4|C/| cosh ^r) |i~3 



RE C {9) 


RE S {9) 




2\U\ 

e {9) + <p * (L ei +L m )(9) 

e^ + ^^ + L^)^), n>l 

Vo (9)+<p*(L m +L eo )(9) 

r ]n (9) + i f *(L Vn _ 1 +L rin+1 )(9), n>l. 



(21) 



These TBA equations have the incidence structure of an infinite chain, where a pair 
of semi-infinite chains of SU(2) invariant factorized scatterings are joined together as 
shown in the following picture: 



£2 ei eo ?7i V2 

■■■ —o o • • o o — ••• 

E c E s 

While the TBA equations are non-linearly coupled equations and are hard to solve 
explicitly, it is rather easy for the cases of XJV(R — > 0) and IR (R —> 00) limits. However, 
the analysis of these two limits depends on the energy-momentum dispersion relations. 

(i) Contrary to the relativistic case where the central charge can still get nontrivial 
contributions in the UV limit because E{6) can be large enough to overcome vanishing 
R in Eq.(13) at large rapidity, the central charge becomes zero for the Hubbard model 
as R — > since P' c {0) and P' s {9) are bound over all values of rapidities. 

(ii) The situation changes in the IR limit. When U is finite, the spinons become 
massive and do not contribute to the central charge because the pseudo-energy 770 (0) ~ 
RE S (6) — > 00 for all values of 0. After decoupling spinon sector, we have only semi- 
infinite chain of holon sector. For a finite value of 0, E c (9) is non-zero and from Eq.(21), 
e o (0) becomes infinite. Therefore, only non- vanishing contribution to the central charge 
comes from — > 00 limit. In this limit, taking a derivative on Eq.(14) and substituting 
P' c with E' c /v c into Eq.(13), we can now invoke the standard TBA analysis to evaluate 
the central charge in terms of the pseudo-energies at = and 00. Note that the Fermi 
velocity in the denominator of Eq.(9) is canceled by the v c in the above substitution. 
The final result for the central charges is 



6 ~ 



n=0 



,-yjOO \ f , y» 



where C(x) is the Rogers dilogarithmic function; 



(22) 



ln(l-*) lnt 



(23) 



t (1-0. 

and we have defined = exp[— e n (oo)] and x° n = exp[— e n (0)]. The TBA equations, 
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Eq.(21), can be rewritten as a set of algebraic equations for x n 's; 



x. 



,00 



~~ 



x 



n 



,00 



(l + :d)*(l + :Ci)*. n^ 1 . 
(1 + 4^)^(1 + n>l. 



(24) 



.0 



ri 



These have solutions x™ — (n + 2) 2 — 1 and 2° = (n + l) 2 — 1. Since x° = o^x, only 

2 

a:^ survives to give c e g = 1 after using £(1) = 

As we claimed in the beginning, the central charge we computed comes from the 
holon sector while the spinon sector decouples. This is what happens for finite U but 
seems valid even vanishing U as far as U > 1/R. In the literature, a discontinuity 
between U — > and U = has been predicted such that if U = the central charge 
will be 2 because the model is nothing but a theory with four free fermions. In our 
analysis, we use the holon and spinon S-matrices which are valid for non-vanishing U. 
Within this validity, our result is what has been noticed in the literature from different 
computations. With this confirmation for the EBler-Korepin S'-matrices, we established 
the S'-matrix program for the non-relativistic models. The most noticeable difference of 
the non-relativistic theories from the relativistic ones is that the IR limit of the former 
models corresponds to the UV limit of the latter. This is the limit where massless 
degree of freedom survives and gives correct central charges. We hope our result can 
be a starting point for the application of the S'-matrices to various lattice problems. 
What we are particularly interested in is to compute correlation functions using the 
form-factor approach |TJJ]. In this scenario, correlation functions of any local operator 
are expressed in terms of the form-factors which can be computed from the exact S'- 
matrices. Although one needs to sum an infinite terms, this can be realized as the sum 
converges very fast. We hope to report this result elsewhere. 
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